Package of facts and theorems for efficiently generating entanglement criteria for 

many qubits 
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We present a package of mathematical theorems, which allows us to construct multipartite en- 
tanglement criteria. Importantly, establishing bounds for certain classes of entanglement does not 
take an optimization over continuous sets of states. These bonds are found from the properties 
of commutativity graphs of operators used in the criterion. We present two examples of criteria 
constructed according to our method. One of them detects genuine 5-qubit entanglement without 
ever referring to correlations between all five qubits. 



For pure states, the notion of entanglement is relatively 
simple to define, even for multipartite cases. A composite 
system is said to be entangled if subsystems are in cer- 
tain states only in reference to each other. For genuine 
multipartite entanglement we need to define properties 
of all but one subsystems via projective measurements 
to have the last one in a pure state. In such situations, 
all reduced density matrices are mixed. 

The problem becomes more involved when we extend 
our consideration to mixed states. We can no longer 
rely on the contrast in purities of the global and reduced 
states. Instead, we need to refer to the operational defini- 
tion. An iV-particle mixed state is genuinely multipartite 
entangled if we cannot obtain it by local actions, classical 
communications, and distributing states that are separa- 
ble with respect to any division of subsystems. 

The difficulties we face in discussing multipartite en- 
tanglement have more than one origin. One problem is 
the relatively large complexity of the question. If we 
want to claim that a state is multipartite-entangled, not 
only do we need to demonstrate that it is not separable 
with respect to any cut, but we also need to show that a 
state is not a mixture of separable states with different 
cuts. The other problem is that multipartite entangle- 
ment can come in various classes of equivalence. For 
example, Diir et al. [l| have demonstrated that a collec- 
tion of 3-qubit W states cannot be transformed by local 
operations and classical communication into a collection 
of 3-qubit GreenBerger-Horne-Zeilinger (GHZ) states, or 
vice versa. The difference between such classes of multi- 
partite entanglement is yet to be fully understood. The 
first step toward such an understanding is to deliver a 
quantifier of entanglement vanishing for one such class 
and not for some other [2] . 

Nevertheless, a number of multipartite-entanglement 
criteria have been developed so far. To give a few exam- 
ples, one of the first criteria is based on the observation 
that states with some separability cannot violate Werner- 
Wolf-Zukowski-Brukner (WWWZB) inequalities 041 as 
strongly as genuinely multipartite-entangled states Q. 



More conventional methods, such as entanglement wit- 
nesses @, H| , have also been employed experimentally [§] . 
For systems with continuous variables, various techniques 
to distinguish multipartite-entanglement have been pro- 
posed [10]. Entanglement quantifiers, such as concur- 
rence, have been generalized to multipartite cases as well 
11]. Laskowski and Zukowski gave a simple criterion to 



distinguish genuine iV-particle quantum correlations and 
those involving fewer particles [12| , using generalized Bell 
inequalities based on a single entry of a density matrix. 
More recently, nonlinear multipartite entanglement cri- 



teria were presented in, e.g., [13l. Il4j. in which the argu- 
ments were based on a simple geometrical fact concerning 
the correlation functions that can be determined by local 
experiments. 

Here we present a useful package of mathematical facts 
and theorems, which allows us to analyze entanglement 
criteria for many qubits in terms of the squares of cor- 
relation tensor elements. We will give an example of its 
applications to construct a 5-partite entanglement cri- 
terion exploiting only four-partite correlations. In this 
way, we can demonstrate "glocal" aspects of some states 
detected by this criterion. On the one hand, entangle- 
ment is globally distributed over all qubits. On the other 
hand, it manifests in correlations between not all subsys- 
tems (hence local). Moreover, the relation between mul- 
tipartite entanglement (ME) and correlations of a given 
order is very loose. For example, all reduced states of 
generalized GHZ states [l6| are purely classical; we can- 
not conclude the presence of ME from them without a 
promise of the purity of the global state. The Smolin 
State [13], on the other hand, has no lower-order cor- 
relations, while it is not multipartite-entangled. In this 
respect, genuinely multipartite-entangled states without 
iV-qubit correlations have been illustrated in 15]. 

Before we start filling out our toolbox, we need 
to formulate our criterion. We discuss the states 
of N qubits. A convenient representation of 
such states is through the correlation tensor, p = 



'W (X) a ,i>, 
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2x2 unit matrix (later occasionally denoted simply as 

1), Oi = X = ( j Q ) , CT 2 = V = { ■ n ) ,<T 3 = Z = 
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Tensor elements can be obtained from an ex- 



-1 / 

periment, T a \, c ... = trp(o a ® o& ® cr c ...). Out of the set 
of 4^ possible operators a a ® o^ ® a c ... we choose some 
subset a and let s denote its elements. [Hereafter, the 
tensor product sign ((g)) will be omitted.] The separabil- 
ity criteria are written in the form of an inequality: 



Q = ]T(tr(p S )) 2 > B 



sc- 



(1) 



When Eq. (JXJ) holds, the state cannot belong to the sep- 
arability class SC. The bound Bsc is dependent on the 
considered SC and a. The question is, how can we find 
Bsc m an efficient way? 

Let us now present our package, which will allow us to 
conveniently obtain Bsc's- It contains two definitions, 
three facts, and four lemmas. 

Fact 1 Two tensor products of element of the Pauli 
group {l,x,y,z}, either commute or anticommute. 

They anticommute whenever they differ at odd number 
of positions, excluding pairs involving the unit operator. 
For example, xxx commutes with xyz and yzl, while it 
anticommutes with xly. 

Lemma 1 Let O = (Oj., ...,O n ) be a vector of traceless 
Hermitian operators with eigenvalues ±1, which all an- 
ticommute with each other. Then for any ( quantum me- 
chanical) state we have ^2^ = i(Oi) 2 < 1. 

Proof: Despite a proof presented in Ref. [l9j|, we pro- 
vide an independent argument here. First we show that 
V ■ O has eigenvalues ±1, provided that V is a real vector 
of the unit length. Consider 



(V ■ Of = V 2 I = I, 



(2) 



where I stands for the unit matrix of a respective dimen- 
sion, d. All the cross terms in the square of V ■ O vanish 
by anticommutativity, and the operator V ■ O is traceless 
by the tracelessness of 0;'s. Any quantum mechanical 
states of this dimension can be written as 



'■O + sP 



(3) 



P = (Pi, Pm) is a vector of hermitian traceless oper- 
ators, which complete an orthogonal basis of operators. 
We demand that tr(Pj) = tr(OiPj) = tr(PjP k ) = for 
j ^ k. The mean value of V-0 with respect to a state p is 
trp(V-0) = V-r, and thus if we set V\\r, we immediately 
see that Y%=i(Oi) 2 = r* < 1. 



Fact 2 Given hermitian operators Oi,...,0„ with eigen- 
values ±1, which all commute with each other, the max- 
imum o/^" =1 (Oi} 2 is n. 

Naturally, this maximum is obtained by taking the aver- 
age with respect to the common eigenstates of O^s. We 
will also use the following fact: 

Fact 3 Linear function obtains its maximum in one or 
more extreme points of its domain. 

The function we will try to maximize is, indeed, a lin- 
ear function (in fact, basically a sum) of variables ai — 
(O'(si)} 2 . The bounds of the domain will also have a 
linear form, being due either to the spectral properties 
of the operators, < et^ < 1, or to the anticommutation 
relations, say, a\ + < 1. 

Lemma 2 Mixing states cannot make the value of the 
left-hand side (LHS) of Eq. {!]) larger. 

Proof: Q is essentially a quadratic, positive semi- 
definite form, which is familiar from Ref. [13|, 



Q = (T,T)= 7 t> (; < 



(4) 
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The "metric" tensor Gy...fc has only entries or 1, and 
Tij...k are elements of the correlation tensor (see Ref. 0). 
Let us now take a mixture p = pp\ + (1 — p)p%, which 
corresponds to T = pT(l) + (1 - p)T{2). Assume that 
(T(1),T(1)) > (T(2),T(2)). Then because of the posi- 
tivity of (T(l) - T(2),T(1) - T(2)) we have 



(T(l).T(l)) > i((T(l),T(l)) + (T(2),T(2))) 

> (r(i),T(2)). 



(5) 



Note that Q would be a mixture of (T(1),T(1)), 
(T(2),T(2)), and (T(1),T(2)) with respective weights 
p 2 , (1 — p) 2 , and 2p(l — p). Hence Q is always smaller 
than (T(1),T(1)) unless p = 1. 



Finally, we define the cut-anticommutativity and cut- 
commutativity. 

Definition 1 Given two operators, C\ = A\B\ and 
Ci = A2B2, which are products with respect to cut A\B, 
we say that they A\B-anticommmute if they anticommute 
on part A, {A\, A 2 } = or on part B, {Bi,B 2 } = 0. 

Definition 2 Given two operators, C\ — A\B\ and 
C2 = A2B2, which are products with respect to cut A\B, 
we say that they A\B-commute if they commute on part 
A, [Ai,A 2 ] = and on part B, [-81,-82] = 0. 
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For example, operators xx and yy commute, but not 
with respect to the cut. Operators x lxl and xylz com- 
mute with respect to all cuts. 

With these definitions, we now present the last two of 
the Lemmas. 

Lemma 3 Consider a set of hermitian operators with 
eigenvalues ±1. {AiBi}ff 1 such that for all i ^ j AiBi 
(A\B)-anticommutes with AjBj. Additionally, every pair 
Ai,Aj and Bi,Bj either commutes or anticommutes. 
Then for all states separable with respect to cut A\B we 
shall have ^2 i=1 (AiBi) 2 < 1. 

Proof: By Lemma 2 we focus only on product states. 
Note that if a unit matrix appears as one of the opera- 
tors Ai it cannot appear as one of the B^s. Without loss 
of generality, we shall call the side, in which the unit op- 
erator does not appear, B, thus all operators Bi anticom- 

mute. Let us group the operators in terms of commuta- 

(k) 

tivity on the side A, so that all operators A\ in the same 
group fc commute with each other, i.e., [A^ , Aj] = 0. 
Each operator of Bi shall also be labeled with k, depend- 
ing on the group its partner A4 is in. We also have 

{Af\A^} = (k^l), 
{Bf\Bf ] } = 0. (6) 

(k) 

Now, because of the anticommutativity of B\ s within 
each group we have 

i 

where (Ak) 2 = max^A^) 2 , i. e., we choose Ak from the 

/ M — 

set of A\ 's. Since operators Ak anticommute with each 
other, we obtain X^-(^fc) 2 — 1 an< ^ hence the conclusion 
of the lemma is proven. 



Lemma 4 A product state with respect to cut (A\B) 
can be a common eigenstate of hermitian operators 
A\Bi, A n B n (being products with respect to this cut) 
with eigenvalues ±1 only if all these operators (A\B)- 
commute with each other. In such a case X)"=i(^^i) 2 — 
n, with the bound saturated for the common eigenstates. 

Under the restriction that all operators either commute 
or anticommute, the proof follows straightforwardly from 
lemma 3. 

We define a set of operators a containing only TV-fold 
tensor products of Pauli matrices. Two such products 
either commute or anticommute. Let our criterion take 
the form 

Q= ]T(trpO) 2 <B SC . (8) 

oea 



Bsc is a bound dependent on the separability class (SC). 
If this bound is beaten, the state p cannot belong to this 
class. The facts and lemmas above validate the follow- 
ing procedure to establish the maxima of Q obtained by 
states belonging to various SCs. We will assign values ei- 
ther 1 or to terms of Q, keeping in mind the restriction 
imposed by the (cut-)anticommutation relations. The 
maximum of Q for a given cut will be given by the in- 
dependence number of the cut-anticommutativity graph 
for operators belonging to a. The independence num- 
ber is defined as the maximal possible number of vertices 
that are disconnected with each other. It is equal to 
the clique number of the cut-commutativity graph with 
the same partition. The clique number is the number 
of vertices in the largest subgraph (clique), all of whose 
vertices are mutually connected. By definition, commu- 
tativity graphs are a complement of anticommutativity 
graphs with the same cut. 

We proceed with an example of a necessary condition 
for the full separability for three qubits studied in [3j-|5|, 

cr = {xxx, yxx, xyx, yyx, xxy, yxy, xyy, yyy}. (9) 

The anticommutation graph of a takes the form of a 
cube with its all four diagonals. Hence the independence 
number of the graph is 4 — the vertices of a tetrahe- 
dron inscribed into the cube are not directly connected 
with each other. If we take states separable with re- 
spect to one of the qubits, new edges appear on the cut- 
anticommutativity graphs — diagonals of a pair of par- 
allel faces. This reduces the independence number to 2. 
Finally, if we insist on the full separability of the state, all 
face diagonals appear as edges, giving the maximum for 
fully separable states equal to 1. This analysis for more 
particles would allow us to easily reproduce the result of 

E 

Another example we present here is the operator set 
that was introduced in Ref. [l8[ to study the Bell in- 
equalities with lower order correlations. Consider 

a = CP(lxxxz) U CP(lzxxz) U CP(lzxzz), (10) 

where 1 denotes the local unit matrix and CP stands 
for a set of cyclic permutations of the argument. This 
symmetry of a will allow us to segregate all bipartite 
states into only three classes. 
First we are going to show 

Proposition 1 For fully separable states, the maximum 
value of Q is 1. 

The construction of a is presented in details in Ref. [l8| . 
The important feature of th elements of a is that for 
every pair of them there is such a position, at which 
one of the elements has x, while the other has z. Hence 
they all (^4|B|C|D|S)-anticommute; a theorem similar to 
Lemma 2 applies here. By Lemma 2, we focus only on 
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fully separable product states. In order to evaluate Q, 
let us write 

(zllll) 2 <a, (2TIII) 2 < (1 - a), 
(lzlll) 2 < b, (I2TII) 2 < (1 - b), 
(llzll) 2 < c , (llzll) 2 < (1 - c), 




FIG. 1: The (A|BCDi5)-commutation diagram for operators 
from from set (|10p . Due to the cyclic permutation symmetry 
of a, analogous graphs apply to all other (4 — l)-cuts. 



Now, let us proceed to more generic cases, demon- 
strating bounds on Q for various states that may have 
different entanglement structures. We shall label the 
qubits with A,B,C,D,E. Figures 1, 2, 3 present graphs 
of cut-commutativity relations between operators from 
a with respect to cuts (A\BCDE), (AB\CDE), and 
(AC\BDE), respectively. Due to the cyclic permutation 
symmetry of cr, any two-partite cut will be equivalent to 
one of these cases. Only those operators that are rep- 
resented by vertices directly connected can be simulta- 
neously assigned values ±1. Thus the maximal value of 
Q is given by the size of the largest clique in the graph. 
One can easily verify that in Fig. 1 there are cliques of 
three vertices. In the other graphs in Figs. 2 and 3 all 
cliques are simply edges. Hence Q < 3 for all bipartite 
separable states. 



(lllzl) 2 <d, (llla-1) 2 < (1 -d), 
(llllz) 2 < e, (llllx) 2 < (1 - e), (11) 



with < a, b, c, d, e < 1. Then, Q is bounded by 



- c)de + (1 - b)cd{\ - e) 
-d)e 



~ c)de + (1 - b)cd(l - e) 
- d)e 



(12) 

I 




FIG. 2: The (AB|CDi5)-commutation diagram for operators 
from set (|10[l . Due to the cyclic permutation symmetry of 
er, analogous graphs apply to all other (3 — 2)-cuts with two 
neighboring qubits at one side. 

On the other hand, Fig. 4 shows a graph of commuta- 
tivity relations between the operators in a. The biggest 
clique found in this graph is of 5 elements, namely zxxzl, 
xxzlz, xzlzx, zlzxx and lzxxz. Hence Q can attain 5. 

We have thus shown the following. 

Proposition 2 For a given by Eq. U0\) . in a given state 
of five qubits p, Q > 1 implies the presence of entan- 
glement in p. In addition 3 < Q < 5 requires genuine 
five-partite entanglement in p. 

Note that the state p given under Eq. (9) in Ref. [3] 
violates a Bell inequality based on a of Eq. (JTUJ) by a 
factor of 1.806. By the Cauchy- Schwartz inequality, this 



I 

Q < abc(l -d) + bcd{l - e) + (1 - a)cde + a(l - b)de + ab(l - c)e 

+ (1 - a)bc{\ - d) + ab(l - c)(l - e) + a(l - b)(l - d)e + (1 - a)(l 

+ (1 - a)b(l - c)(l - d) + a(l - - c)(l - e) + (1 - a)(l - 6)(1 

+ (1 - a)(l - c)d(l - e) + (1 - b)c(l - d)(l - e) 

< abc(d - 1) + bcd{l - e) + (1 - a)cde + a(l - b)de + ab(l - c)e 

+ (1 - a)bc(l -d) + ab(l - c)(l - e) + a(l - 6)(1 - d)e + (1 - a)(l 

+ (1 - a)b(l - c)(l - d) + a(l - - c)(l - e) + (1 - a)(l - 6)(1 

+ (1 - a)(l - c)d(l - e) + (1 - b)c(l - d)(l - e) 

+ abcde + (1 - a)(l - b)(l - c)(l - d)(l - e) 

= 1. 




5 



i.\n.\ 




FIG. 3: The (AC|BDi5)-commutation diagram for operators 
from set (|10|) . Due to the cyclic permutation symmetry of a, 
similar graphs apply to all other (3 — 2)-cuts with two non- 
neighboring qubits at one side. 



of matching vertices of iV-cube graph into groups of 1, 2, 
A... It is known that the number of such matching grows 
superexponentially with respect to N [20]. Therefore, our 
toolbox is highly efficient in producing a large number of 
multipartite entanglement criteria in a straightforward 
and transparent manner. 

The project was part of the EU project QUASAR. MW 
is supported by the Foundation for Polish Science (HOM- 
ING PLUS Program). KM is grateful for the support of 
the JSPS Kakenhi (C) no. 22540405. 




FIG. 4: The commutation diagram for operators from set 
(Ell- 



implies Q > 3 (see Ref. [18j|). Thus p is genuinely five- 
partite-entangled despite the absence of 5-partite corre- 
lations. 

In conclusion, we have shown a useful toolbox for for- 
mulating (multipartite) entanglement criteria. These cri- 
teria take the form of an inequality, with a separability 
class-dependent upper bound on one side, and an incom- 
plete sum of squared correlation tensor elements on the 
other. Importantly, our method eliminates the necessity 
of finding the optimal state maximizing the expression as 
in Refs. [11 Bj]. 

We have also shown an example with a set of operators 
that was studied in Ref. 18[ in the context of the Bell 
inequalities with lower-order correlations. Interestingly, 
we are able to test TV-partite entanglement without using 
iV-partite correlations. It is important to note that non- 
classical lower-order correlations are characteristic of W- 
type entangled states, whereas they do not appear in 
GHZ-like states. Hence we can additionally argue that 
the entanglement criterion given in example 2 does not 
detect 5-partite GHZ entanglement. 

Lastly, let us stress that we have presented a family of 
multipartite entanglement criteria, rather than a single 
one of them, given by all possible sets of observables cr's. 
Those cr's with only three observables at each side (say 
x, y, and 1) can be obtained by considering the problem 
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